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Abstract 


A spectral analysis of the vertical positions and velocities of 374 open star clusters (OSCs) was carried out. We 
took these OSCs from the Hunt and Reffert catalog; they have an average age of about 10 million years, and are 
located on the galactic plane XY in a narrow zone inclined by 25° to the galactic axis Y. The following estimates of 
the parameters of the Radcliffe wave were obtained: (a) the maximum value in periodic perturbations of vertical 
coordinates Zmax = 92 + 10 pe with the wavelength of these perturbations A,= 4.82 + 0.09 kpc; (b) the 
maximum value of the velocity of vertical disturbances Wnax = 4.36 + 0.12 km s7! with disturbance wavelength 
Aw = 1.78 + 0.02 kpc. Note that the results of the vertical velocity analysis are first-class in accuracy and 


completely new. 
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1. Introduction 


The Galaxy’s disk is far from dynamic equilibrium. 
Perturbations of the coordinates and velocities of its compo- 
nents are detected at various scales. The causes of such 
disturbances are both internal and external factors. Internal 
reasons are responsible, for example, for the presence of a 
large-scale galactic spiral structure and a central bar (Lin & 
Shu 1964; Debattista & Sellwood 2000). External gravitational 
influence is usually sought to explain the effect of large-scale 
curvature of the disk in the vertical direction (for example, 
Bekki 2012). Other disturbances in the spatial velocities of 
stars, in particular vertical ones, have also been discovered, 
which are not associated with the spiral structure (for example, 
Widrow et al. 2012; Antoja et al. 2018). 

In the nearest circumsolar neighborhood with a radius of 
about 2.5 kpc, a Radcliffe wave was recently discovered (Alves 
et al. 2020), the study of which is the main goal of this work. 
This wave was first identified by Alves et al. (2020) when 
studying a large sample of molecular clouds. The main feature 
of the Radcliffe wave is that in the galactic plane XY the clouds 
are elongated “in a string,” which is inclined at an angle 
B ~ 30° to the galactic axis Y, and their vertical coordinates are 
Z that have a wave-like character. The maximum value of the Z 
amplitude is observed in the immediate vicinity of the Sun. 

The behavior of vertical coordinates characteristic of the 
Radcliffe wave is manifested in the distribution of interstellar 
dust (Lallement et al. 2022), molecular clouds (Zucker et al. 
2022), masers and radio stars (Bobylev et al. 2022a, 2022b), T 
Tauri stars (Li & Chen 2022), massive OB stars (Donada & 
Figueras 2021; Thulasidharan et al. 2022), as well as young 
scattered stars clusters (Donada & Figueras 2021). 


The vertical velocities of objects in the Radcliffe wave have 
not yet been studied well enough. Donada & Figueras (2021) 
apparently found the first connection between vertical coordi- 
nates and vertical velocities from an analysis of 11 open star 
clusters (OSCs). Thulasidharan et al. (2022) showed that the 
maximum amplitude of vertical velocities of OB stars in the 
Radcliffe wave can reach 3-4 km s™!. Bobylev et al. (2022a) 
found a relationship between the vertical coordinates and 
vertical velocities of masers and radio stars. 

A number of hypotheses for the origin of the Radcliffe wave 
have been proposed. For example, Fleck (2020) links the origin 
of the Radcliffe wave to the Kelvin-Helmholtz instability. 
Marchal & Martin (2023) connect the origin of the Radcliffe 
wave with the evolution of the North Celestial Pole Loop 
(NCPL). In their opinion, a narrow chain of material (from 
which the Radcliffe wave was formed) arose at the outer 
boundary of the NCPL as a result of the impact of shock waves 
from several supernovae and their stellar winds. A hypothesis 
has also been put forward about the impact of an external 
impactor on the galactic disk (Thulasidharan et al. 2022; 
Oseguera et al. 2023). Such an impactor could be a dwarf 
galaxy or a globular cluster. The passage of such a massive 
body through the galactic disk can cause wave disturbances 
propagating across the disk. Each of the listed hypotheses has 
advantages and disadvantages, but there is currently no 
generally accepted hypothesis for the formation of the 
Radcliffe wave. 

For a detailed study of the distribution and kinematics of stars, 
measurements of such characteristics as trigonometric parallaxes, 
proper motions, and radial velocities of stars are important. 
Recently, such measurements have emerged as a result of the 
Gaia space project (Gaia Collaboration et al. 2016). In the latest 
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Figure 1. Distribution of selected OSCs younger than 25 million years in 
projection on the galactic plane XY; the selection zone is shown by shading; the 
Sun is located at the point with coordinates (X, Y) = (8.1, 0) kpc. Roman 
numerals indicate the position of the spiral arms: I—a segment of the Perseus 
arm; II—a segment of the Carina-Sagittarius arm. 


published version of the project, in the Gaia Data Release 3 
(DR3) catalog (Gaia Collaboration et al. 2022), the proper 
motions of about half of the stars in the catalog were measured 
with a relative error of less than 10%, and their radial velocities 
were measured for a large number of stars. 

The presence of three projections of star velocities allows us 
to perform a full three-dimensional analysis of their kinematics. 
This is especially important when studying the Radcliffe wave, 
since only the spatial distribution of its components has been 
well studied. All the more interesting are the kinematic data not 
about single stars, but about OSCs, where more accurate 
average values have been calculated for a large number of 
cluster members. One of the most extensive modern kinematic 
databases on galactic OSCs is the Hunt & Reffert (2023) 
catalog, which contains the parameters of 7200 OSCs identified 
according to the Gaia DR3 catalog. In this work, we apply 
spectral analysis to the coordinates and vertical velocities of 
OSCs to clarify the features of the Radcliffe wave using a large 
sample of young OSCs from Hunt & Reffert (2023). 


2. Method 


Two coordinate systems are used: (1) a rectangular 
heliocentric coordinate system, in which the x axis is directed 
from the Sun to the center of the Galaxy, the y axis coincides 
with the direction of rotation of the Galaxy and the z axis is 
directed to the north galactic pole; (2) galactocentric rectan- 
gular coordinate system X, Y, Z, in which the direction of the X 
axis is opposite to the direction of the x axis, and the directions 
of the Y and y axes, as well as the Z and z axes, coincide 
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between themselves. So, in these two coordinate systems, only 
the directions of the x and X axes differ. 

The orientation of the Radcliffe wave relative to the y and Y 
axes differs only in sign. In the coordinate system xy, let us 
make the transition to the hatched axis y’ by rotating the system 
by an angle @. As a result, we obtain a connection between the 
coordinate y’ and the coordinates x, y in the form: y'= 
ycos 8 + xsin p. 

To study the periodic structure in the coordinates and 
velocities of stars, Bobylev et al. (2022a, 2022b) proposed to 
use spectral analysis based on the standard Fourier transform of 
the original sequence z(y’) : 


FON = f z(y) edy! 
= U(A) + jV (X) = AA) ef, (1) 


where A(\) = U?(A) + V7(A) is spectrum amplitude, and 
(àA) = arctan(V(A)/U(A)) is spectrum phase. We limit the 
resulting spectrum to the wavelength range that most accurately 
describes the original data. Essentially, it coincides with the 
main peak (lobe) of the calculated spectrum in the wavelength 
range from Amin tO Amax (within these boundaries the spectrum 
smoothly decreases starting from the maximum value, and 
outside —starts to increase). 

The desired curve, which approximates the original data, is 
calculated using the inverse Fourier transform formula in the 
wavelength range we have determined: 


Try! 


Amax 
z(y) = 2k f . AQyeos( 27 


+ po Jaa, (2) 


where the coefficient k is found from the condition of minimum 
residual. 


3. Data 


The catalog of Hunt & Reffert (2023) includes 4780 OSCs 
already known before, and 2420 of the total number found are 
new candidates. For all OSCs in the catalog, ages, lifetimes, 
distances, and kinematic characteristics were estimated. An 
important advantage of the Hunt & Reffert (2023) catalog is the 
availability of average radial velocities for a large number of 
OSCs, many of which are new. Moreover, the average radial 
velocities of all OSCs were calculated exclusively from the 
data of the Gaia DR3 catalog. 

In this work, clusters were selected in a zone limited by two 
lines on the xy plane: 


x = ytan 8 — 0.65, 
x = ytanG + 0.10, (3) 


where the angle G=25° was selected so that the wave 
manifests itself in the best possible way. The center of the 
selection zone does not pass through the origin of coordinates, 
but cuts the x axis at the point x,; therefore, the transition to the 
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Figure 2. OSC coordinates z as a function of distance y’ (a) and their power spectrum (b); the periodic curve shown by the thick line reflects the results of spectral 


analysis. 
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Figure 3. Vertical velocities of OSCs W as a function of distance y’ (a) and their power spectrum (b); the periodic curve shown by the thick line reflects the results of 


spectral analysis. 


y’ axis was performed according to the formula 
y! = ycos3 + (x — x)sin 8, (4) 


where xı = — 0.28 kpc. OSCs younger than 25 million years 
were taken. The selection zone contained 374 OSCs with an 
average age of 9.7 Myr. The distribution of these OSCs in 
projection on the galactic plane XY is given in Figure 1. 


4. Results and Discussion 


The positions of the young OSCs we selected in the galactic 
plane XY were projected onto the y’. axis. Next, a spectral 
analysis of the vertical positions and vertical velocities of the 
OSCs was carried out in this coordinate system. Based on the 
distribution of open space in the vertical direction, the 
following estimates of Zmax and A were obtained: 


Zmax = 92 + 10 pe, 
A = 4.82 + 0.09 kpc, (5) 


and the value Zmax is achieved at y’ = —0.28 kpc. Based on the 
vertical velocities of the OSCs, an estimate of the maximum 
value of their disturbance velocity Wmax with the wavelength of 
these disturbances À was obtained: 


Woax = 4.36 + 0.12 km s7}, 
A = 1.78 + 0.02 kpc, (6) 


and Wmax is achieved at y’ = —0.71 kpc. The results of the 
spectral analysis are shown in Figures 2 and 3. 

The value of the wavelength A found in solution (5) differs 
markedly from the results of other authors (where A ~ 2.5 kpc 
is usually the case), but our estimates of the W velocities 
(solution 6) are first-class in accuracy and completely new. 

In the work of Bobylev et al. (2022a), two samples of young 
stars, namely a sample of masers with measured trigonometric 
parallaxes and a sample of T Tauri stars were studied. In both 
cases, these are protostars and very young stars of various 
masses. After appropriate selection of stars in a narrow band, it 
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was shown that the Radcliffe wave manifests itself in the 
positions and velocities of both masers and T Tauri stars. Based 
on the spectral analysis of the masers, the following estimates 
were obtained: Zmax = 87 +4 pe and A=2.8+0.1 kpc, 
Wax = 5.1 + 0.7 km s | and the wavelength found from the 
vertical velocities À = 3.9 + 1.6 kpc. The following estimates 
were obtained for T Tauri stars: Zmax = 118 +3 pe and 
A=2.0+0.1 kpc. 

Analysis of various populations of young objects belonging 
to the Radcliffe wave does not give an unambiguous picture of 
the behavior of vertical velocities along the wave. Tu et al. 
(2022) found no significant changes in the vertical velocities of 
young stars depending on the position along the wave. These 
authors concluded that such changes could be less than a few 
kilometers per second. Li & Chen (2022), studying the proper 
motions of young stars that have not reached the main sequence 
stage, found a phase difference of about 27/3 between the 
waves of vertical coordinates and velocities. This conclusion 
contradicts the results of the analysis of W velocities obtained 
by Donada & Figueras (2021), Thulasidharan et al. (2022), 
Bobylev et al. (2022a, 2022b). This conclusion also contradicts 
the results of this work. Indeed, based on the position of the 
main minimum in Figures 2 and 3, we cannot talk about a 
significant shift in the phases of the two waves. 


5. Conclusion 


The study of the spatial and kinematic characteristics of the 
Radcliffe wave was carried out using a sample of young OSCs. 
For this purpose, we selected clusters younger than 25 Myr 
from the modern catalog of Hunt & Reffert (2023). The 
selection of OSCs was carried out from the near-solar 
neighborhood with a radius of about 3 kpc, in a narrow strip 
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located at an angle of 25° to the galactic axis Y. There were 374 
such clusters, which have an average age of 9.7 Myr. It is 
shown that there is a periodicity characteristic of the Radcliffe 
wave, both in the vertical coordinates and in the vertical 
velocities of the OSCs. 

To evaluate the characteristics of the Radcliffe wave, we 
used spectral analysis. As a result, the following estimates were 
obtained: (a) maximum value of the coordinate z, Zmax = 92 + 
10 pe and disturbance wavelength A= 4.82 +0.09 kpc; (b) 
maximum value of vertical velocity disturbances Wmax = 
4.36 + 0.12 km s_' and the wavelength of these disturbances 
A= 1.78 + 0.02 kpc. 
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